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Abstract 



Green function and linear integrals of motion for a charged particle moving in elèctric field are 
discussed. Wigner function and tomogram of the ststionary states of the charge are obtained. Con- 
nection of quantum propagators for Schròdinger evolution equation, Moyal evolution equation and 
evolution equation in tomographic probability representation for charge moving in elèctric field is 
discussed. 



Introduction 



o 
o 

o 

C^he quantum mechanics formalism is based on the notions of wave function |J ^ [| |], [5] and density 
matrix |J. The possible interpretations of quantum properties including the interpretation in terms of 
^Jiidden variables have been suggested in JF], |^| and analysis of this approach is presented in [TO, O. The 
Qomographic map [E| gives the possibility to formulate quantum mechanics using the notion of probability 
^nd, in principle, avoiding the notion of wave function fï^, [ij], [ï^, |Tj| . The attempts to construct quantum 
O-formalism similar to classical one done in the spirit of Moyal approach 
Ajvere made using quasidistributions of different kind [I 



IQI |2T], p2| . The evolution equation of such 
O^uasidistributions |23], |Ï3| provides the complete information on behaviour of the quantum systems which 
np f - equivalent to the information contained in von Neumann equation for density matrix ||. There are 
Q^veral modern quantum problems |23|, |25|, |27], |28|, ^9], |3Ü], |3Ï|, [32], |33], [34]] related to trapped ions, 



based on Wigner function [18 



^hotons in cavity, quantum computers as well as traditional ones |35], |36| which can be considered in 
OT"amework of the probability representation of quantum mechanics. In this work we analyse one of these 
^í'oblems for charge moving in elèctric field. 

L| ■ One of the solvable problems of quantum motion is the problem of charge moving in elèctric field. It is 
• Hjmportant problem and we will discuss this problem in framework of tomographic probability representation 
çàf quantum states [g?], gj| H M> S Él 0> @ ■ 



,_ ctf \ The probability representation of quantum mechanics was introduced in [|T3], [Ï4| . The Wigner function 



18] of quantum states is the informative characteristic of the state. The optical tomography scheme |24 



was used to reconstruct the Wigner function of photons. Extension of the optical tomography method 



containing two real parameters was suggested in |L2| . The tomographic probability distribution of quantum 



state was employed in [TJj] to replace the wave function and density matrix by the probability density 
which determines the state. This approach to description of quantum states created new formalism of 
quantum mechanics similar to classical statistics. The evolution |Ï3 , 14 , |43l , [Ï6| and energy levels JÏ5, 
[Ï6| in tomographic representation correspond to Moyal formulation of quantum mechanics fTTf . The 
probability representation was introduced for spin in [[Ï4], ^5], £Ï7|, |l| . The Pauli equation in tomographic 



representation was obtained in |49| . There are quasidistribution functions which correspond to density 
matrix in different representations like Wigner function |1§, Glauber-Sudarshan P-function JÏ9| and p0| , 
Husimi Q-function [^TJ. The tomograms have specific property of Standard probability distribution and 
they completely describe the quantum states. Many physical systems are described by the quadratic 
Hamiltonians pÜf . The charge moving in homogeneous elèctric field is important example of such systems. 



1 



The Gaussian and Gauss-Hermite states for such system were studied using time-dependent integrals of 
motion linear in position and momentum in ^TJ. The wave functions of Gaussian states and number states 
can be obtained for the charge moving in elèctric field in explicit form. The tomographic description of the 
states is the main goal of our work. We will obtain tomograms and propagator for the evolution equation 
in the probability representat ion for the charge moving in elèctric field. We will also discuss the connection 
of Green function for Schròdinfer equation and the propagator. 



2 Green function and integrals of motion 

First we discuss the Green function of Schròdinger evolution equation for the wave function in position 
representation 

itj ^A = H(tMx,t). (1) 

Unitary operator U(t) is called evolution operator of this equation if for any solution ip(x,t) expressed in 
terms of the state vector \ip,t) one has the relation 

hM) = £(í)|V,0). (2) 

From this equation follows that at the time moment t = the evolution operator coincides with the 
identity operator: 

Ü(0) = T. (3) 

The matrix element of the operator U(t) in position representation is called "Green function" of the 
Schròdinger evolution equation (J^) 

G(x,x',t) = (x\Ü(t)\x'}. (4) 
In position representation the equation (|J) can be rewritten in the integral form 

i{>(x,t) = / G{x,x',t)ifj{x',0)dx'. (5) 



This equality means that the Green function is the kernel of the evolution operator. It is obvious that for 
the initial time moment the Green function coincides with Dirac delta-function, i.e., 

G(x,x',0) = 5(x-x'). (6) 

It is easy to obtain the evolution equation for the Green function G(x,x',t). Substituting expression (|) 
into Schròdinger equation for the wave function we get 

ih^(x,t) = ih J ^(x,x',t)ip(x',0)dx' = H(t) [x] i/)(x,t) = J H(t) [x] G(x,x',t)ip(x / ,0)dx\ (7) 

where notation H\ x i means that the energy-operator H acts on the variable x of the Green function. Since 
the above equations are vàlid for arbitrary wave function we obtain the following equation for the Green 
function 

dG ~ 

ih— (x, x', t) = H(t) [x] G(x, x', t). (8) 
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Now we discuss the notion of quantum integrals of motion. The operator I(t) is called the integral of 
motion if its mean does not depend on time, i.e. 

(T(t))^(^t\T(tM,t) (9) 

for any solution \ip,t) of the Schròdinger equation. It means that 

ü %)) =^,t\Tm,t) = o. (io) 



dt\ I v dt 

Since the total time derivative can be expressed in the form containing the term with the commutator 

the operator I(t) is integral of motion iff it satisfies the evolution equation of the form 

ift^ = [£,/(*)], (12) 

in other words, the operator I(t) is the integral of motion iff its total time derivative is equal to zero: 

dT(t)_dT(t) i m? 

By definition of the evolution operator U (t) one has the equalities 

foM) = Ü(tM0), = (^0\Ü + (t) = {^0\Ü-\t). (14) 

We used the property of unitarity of the evolution operator which takes place for Hermitian Hamiltonians. 
By definition of the integral of motion I(t) we have also the equalities 

(^0|/(0)|^,0) = (^í|/(í)|^,í) = ty,0\Ü-\t)T(t)Ü(t)\il>,0), (15) 

from which one gets 

7(0) = U- l {t)í{t)U{t) or T(t) = ^(í)f(0)^ _1 (í). (16) 
This mean that the evolution operator and the integral of motion satisfy the operator equation 

T(t)ü(t) = ü(t)T(o). (i?) 

The matrix elements of the equal operators must be equal, i.e., 

%)Ü(tj) =(Ü(t)T(0)) , (18) 



The above equation is written in position representation. Note that for any operator / and function ip(x) 
the following equality takes place 

(fi/>)(x) = J f(x,x')i>(x')dx', (19) 



which is an analog of finite-dimensional equality 

n 

/(Oi = (20) 

where / : L — > L is linear operator in n-dimensional linear space L, l G L is n-vector, and Zj are matrix 
elements of operator / and components of vector l, respectively, with respect to some fixed basis of linear 
space L. Equation (|i8|) can be presented in the form 



I(x,y,t)G(y,x',t)dy = I(t) [x] G(x,x',t) = / G(x,y,t)I(y,x',0) dy = / I(y,x',0)G(x,y,t) dy 

J J (21) 

P(y, x\ 0)G(x, y, t) dy = ?(0) M G(i, x', t). 



Here P implies the transposed operator. Thus, we obtain system of two equations for Green function in 
case of physical system with one degree of freedom 

T(t) [x] G(x,x',t) = í i (0) [x , ] G(x,x',t), (22) 

where the operator I(t) has two components which are independent integrals of motion. This system of 
equations together with equation (ffi) and the initial condition (||) provides the complete system of equations 
determining the Green function. The normalization constant for the Green function can be found using 
the nonlinear equation for the evolution operator 

Ü(t) = Ü(t-T)Ü(r). (23) 

This operator equality gives the nonlinear integral equation for the Green function 

G{x,x',t)= \ G{x,y,t — r)G(y,x' \r) dy. (24) 



This equality is equivalent to the initial condition @ for the Green function. 

The considered general properties of the quantum integrals of motion and Green function of Schrodinger 
evolution equation will be studied for partial quantum system. We consider a particle in a uniform elèctric 
filed F. The Hamiltonian H of this particle reads 

- p 2 

H = — Fx. (25) 

2m v ' 

If it were the classical particle we would have the trajectory of the particle in the phase space of the form 

p Ft 2 

x = x H 1+ — , (26) 

m 2m 

p = Po + Ft (27) 

here x is the coordinate and p is the momentum of the classical particle, respectively. Hence, the inverse 
expressions 

P Ft 2 

XQ = x -^t + —, 28 
m 2m 

p =p-Ft (29) 



4 



provide the integrals of (classical) motion. It is easy to check that operators 

~ ~ P Ft 2 

x Q = x 1+ - — , (30) 

m Ira 

p =p-Ft (31) 

are integrals of (quantum) motion, i.e. they satisfy the equation (|l2f) . These operators are obtained by 
means of quantization procedure applied to the classical integrals of motion. The physical meaning of the 
discussed integrals of motion is explained by the fact that these operators describe the initial point in the 
particle phase space. 

For Green function G(x,x',t) we have the system of two equations j5(| 

x (t) [x] G(x,x' ,t) = x'G(x,x',t), 

„t 1 \ * dG , / \ ( 32 ) 
p Q (t) [x] G{x, x,t) = ™Qg\*, x i *)> 

or for the charge moving in elèctric field 



%fxt d Ft i ! i , 

x H — — I G(x, x , t) = x G(x, x , t), 

m ox 2m ' 



(33) 



(d \ dG 

ih— — h Ft 1 G(x, x' t) = —ih— — (x, x' t) 
ox ) ox' 

with the initial condition @ 

G(x,x',0) = 5(x-x'). (34) 

The system of equations is obtained using the quantum integrals of motion. To solve this system we try 
to find the Green function in the Gaussian form, i.e. 

G(x, x', t) = N(x\ t)e a{x '' t)x2+b{x '' t)x , (35) 

where N(x',t), a(x',t), and b(x',t) are unknown functions. Substituting this expression into the first 
equation of the system (p3|) we have the equality 

x + — Í2a(x', t)x + b(x', t)) + — = x', (36) 
m V / 2m 

from which we obtain the functions a(x',t) and b(x',t) in the form 

, , . im , . , , ira . iFt , . 

«t) = -, b ^, t) = -- x ' + - (37) 



To find out the function N(x',t) we substitute ( pop into the second equation of the system ([$3|). As a 
result we obtain the equation for this function: 

dN i fra . 1 \ 

— — = — ( — x H — Ft ) N. 38 

dx' h V t 2 J K J 

To solve this equation we try to find the function N(x',t) in the Gaussian form, i.e., 

N(x', t) = n(t)e a{t)x ' 2+mx ' . (39) 
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Substituting this expression into previous equation we obtain the functions a(t) and (3{t): 

im „. . iFt , . 

°w = ss- " (t) = ÜT ' (40) 



Thus the explícit form of the Green function G(x, x', t) reads 

{ 



G(x,x',t) =n(í)exp<j ~ x 'f + l^( x + x ') ^ ( Jl > 



where the preexponential factor is unknown function of time. To find out the function n(t) we substitute 
this expression in the equation (|$D and get the differential equation of the form 

dn ( 1 iF 2 t 2 \ . , 

+ -5ZT- K (42) 



The solution to this equation reads 

This solution contains unknown constant c and for the Green function we have the explicit expression 

G (x,x', t ) = 7! exp|-(x-,0 + _( I + a 0-_j. (44) 

The constant c can be obtained from the nonlinear equation (0). If we substitute the expression (|44]) into 
this equation we obtain the equality 



2nih , , 

l = c\/ ■ 45 

m 



The above equality determines the normalization constant of the Green function 



m 



Finally we have the following expression for the Green function 



(46) 



/ m I im , , 2 iFt .. iF 2 t 3 1 

G(x,x,t) = J— — exp<^ — (x-x) + —( x + x)-—-r}. (47) 



)< — — (X — X) H — (X + X ) ; >. 

\2hV ' 2h y ' 24mhj 



which has the Gaussian form. The expression in curle brackets is the classical action S(x,x',t) multiplied 
by the factor i/H. For F = we get the Green function for free motion. 

The equation for density operator p(í) is similar to Schròdinger eqution for the wave function 

ihp(t) = [H(t),p(t)}- (48) 
From this equation it is easy to obtain equation for matrix elements of density matrix p(x, x', t) = (x\p(t) \x') 

^(x, x\ t) + l - (l·l [x] p(x, x', t) - H\ x , ]P {x, x\ i)) = 0. (49) 
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The Hamiltonian of moving particle with potential V(x) reads 



S = £ + V W , (50) 



and equation takes the form 



!<*■ 4) " " J*) " s K> ~ V{X,) Y) «*• •'· «> = (51) 

The equation for density matrix admits the solutions in factorised form if)(x, t)ip*(x' , t) which describe the 
pure states of the quantum system. 

The Schròdinger equation for stationary state of the particle in the elèctric field reads 

^ + -( £ + Fx)i> = . (52) 

It is well known that for all energy vàlues £ there exists the solution ip£Í x ) °f equation (|52| ) given by a 
formula 



where is the Airy function: 



oo 

3 \ i r ( /„,3 



= —j= i cos ( — + ux J rfw = _ / exp |ï ( — + ) /> du, (54) 





and the coefficient in front of the function $ is necessary for normalization of ip£ to 5-function: 

il>e(x)i/i e (x) dx = 5 (£' - £) . (55) 

Plot of the Airy function is presented in Figure (JH). Let us introduce notations 

(2m) 1/3 f2mF\ 1/3 £ . . 

^= ^1/6^/3 ' a = -( v -^J ' £ = a F- (56) 

In this notations ips reads ips(x) = (ax + e). 



3 Weyl representation 

The density operator of quantum state can be considered in different reprresentations. One of the repre- 
sentations is Weyl representation. In this representation each operator is described by its symbol. The 
Weyl symbol of the density operator is called Wigner function. 

Now we obtain the equation for the Wigner function. The Wigner function W(q,p) is expressed in 
terms of matrix elements of density operator as [jj] 

W(q,p)= í p(q + ^ q -^ e -^ u du. (57) 
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Figure 1: Airy function 



To obtain the equation for the Wigner function we have to know the correspondence between operators 
acting on the functions of variables x, x' and the same operators acting on the functions of variables q, p. 
The operators /[ XjX '] and F[ gjP ] are called corresponding to each other 

f[ x ,x'] < — ► F{q, P ] (58) 

if the following equivalent relations are vàlid for any corresponding to each other density matrix p(x, x') 
and the Wigner function W(q,p) of the same quantum state 



f (W) (ç + l^-f) e~ï pu du= (F [q , p] w) (q,p), 



(59) 



and 



p(x,x') = ^jw (^,p) et«*-Vdp=ï 
(w) M = \f (F M W) ei^'Up. 
Now we obtain operators corresponding to the operators x, x', d/dx, and d/dx'. 

J (xp) (q+l,q-l)e-^du = f + e-ïvu du 



(60) 



(61) 



It means that multiplication of density matrix by the coordinate x provides the action on the Wigner 
function of the form 

xp(x,x')< ► í q + —— \ W(q,p). (62) 
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Analogously one can obtain that 



x'p(x,x')< ► ( q- —— ) W(q,p). (63) 



To obtain operators acting on the Wigner function corresponding to operators x and x' we used the 
equation fl59|) . Now to obtain operators acting on the Wigner function corresponding to operators d/dx 
and d/dx' acting on the density matrix we use the equation (|60D and get 



9p ^x') = XJ(~ + 1 JP) W (^,p) et*-*Up. (64) 



dx 2nh 
By definition of corresponding operators we have 



d , (\ d i 

—p(x,x') < — ► -— + -p ) W(q,p), (65) 



dx \ 2 dq h 

and analogously 

é P{X > X,) "{líq-l P ) W{q > P) - (66) 

Now using the obtained correspondence rules and the von Neumenn equation (|5Ï| ) for the density matrix 
we can write the equation for the Wigner function, which was obtained by Moyal |Ï7 



8 \ P ,.) + ^( í ,P,<) + iM í+ ?^-^ í 4^]l%M) = 0. (67) 



dt m dq h\ \ 2 dp J \ 2 dp 

In the case of charged particle moving in elèctric field the potential energy V(x) = —Fx, and one gets the 
Moyal equation in the form 

• p dW 1 dW , . 

m dq 2 dp 

The Wigner function W(q,p,t) can be found as solution to the Moyal equation. It has the form 

W(q,p,t) = W (x ,p ) (69) 

where Wo(q,p) is the initial Wigner function and xo and po are given by equations (|28|) and (p9[), respec- 
tively We get 

W(q, p,t)=W (q--t + ^,p-Ft). (70) 
\ m 2m J 

For Wigner function one can introduce the propagator T(q,p, q',p', t) of Moyal equation 

W(q,p,t) = J r(q,p,q l ,p l ,t)W(q',p , ,0)dq'dp'. (71) 

In the case of the charge moving in elèctric field we have explicit expression for the propagator of Moyal 
equation 

( P F t 2 \ 

r(q,p,q',p',t) = 5(x -q')S(p - p') =S[q 1 - q + — )6 (p - p - Ft) . (72) 

V m 2m I 
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Let us consider example of initial state which is described by the wave function 



.meu meu ., . 

Mx) = \J^h exp \-^h x >■ (73) 

This state can be specially prepared by means of a spring which is removed at the time moment t — 0. 
The corresponding Wigner function has the initial form 



rau 



P 2 \ 



^ (ï ,p) = 2exp|--^-^|. (74) 
The evolution of the Wigner function is given by formula 



m·jj 2ut 1 ( 2 1 \ L > Ft 3 üü 



W{q,p,t) = 2exp^ — q + —qp- — [ut + -)p - — p 

n n hm \ u J hm 

l·'t'.c FH 2 fujt 2 1 \ 2/7 



-r + - + 



h hm \ 4 w / hmtu 



(75) 



Now we calculate the Wigner's function of the stationary state of the charged particle with given energy 
3|). By definition, Wigner function of the state is related to the wave function ips(x) by the formula 



W e (q,p) = J^s (q + |) Pe (q ~ f ) e^du. (76) 
Using expressions ( |53D and (|54]) it is easy to show that 

1 ,Hrn~^ í ( 1 \ 2/3 9 ímF\ 1/3 / m \V3 



(77) 



Plot of the Wigner function is shown in Figure @. 

4 Tomographic representation 

The tomogram of the quantum state is expressed in terms of Wigner function (h = 1) ||12|| . 

wiX^u) = í W( g ,p)e^-^ (78) 
J (2tt) 

The inverse relation determines the Wigner function in terms of the tomogram 

W(q,p) = f w(X,»,u)e^^ (79) 

For pure state described by wave function ip(x) for tomogram one can obtain the formula in terms of this 
wave function 



w(X, /i, v) 



2tï\v\ 



r(./y)<>xp^ |V - l —y \ d V 
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15 

Figure 2: Wigner function 
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To write the equation for tomogram [14] we have to know the correspondence between operators acting 
on Wigner functions depending on the variables q, p and the same operators acting on the tomograms 
depending on the variables X, fi, v. Operators -F[ 9lP ] and G[x,n,u] are called corresponding to each other 



[q,p] *- 



G 



51) 



if the following equivalent relations are vàlid for any corresponding to each other Wigner function W(q,p) 
and tomogram w(X, fi, v) 

dk dq dp 



w{X,fi,v)= W{q,p)e- lk{x - m - up) 



F[ q>P ]W) (q,p)e 



and 



W{q,p) = / w(X,fi,u)e 



(27T) 

ik(x-nq-vp) dk dq dp _ 

(27T) 2 

Hx-w-vp) dXdfidv 



(82) 



G[x,n, v ]w) (X,fi, v) 



2tt 



F [q , P ]W) (Q,P) = / {G[x,,M w ) {X,fi,v)e^ x -^ 



dX dfi du 
2^ ' 



53) 



Using the same method which we used for finding the corresponding operators in position and Weyl 
representation we can get the following relations 



qW(q,p) 
pW(q,p) 







dX ) 


dfi 


d > 




dX J 





w(X,fi, v) 
w(X, fi, v) 



d_ 
dq 
d_ 

dp 



W{q,p) 



d 
d 

W(q,p) < > v—w{X,fi, v). 



54) 
55) 



The action of the operators onto homogeneous tomograms must give the homogeneous function. This 
condition gives the possibility to avoid the ambiguity related with the fact that Wigner function depends 
on two variables and the tomogram depends on the three variables. 



The evolution equation for tomogram in which we reconstruct the Plank constant reads |14 

-i 



dw dw i 

~dt~^~dv~ + h 



V\- 



d 

dX 



1 d ihv d 
~dfi + TdX 



V 



d 

dX 



d_ 
dfi 



ihv d 



w = 0. 



56) 



2 dX 

In our case of charged particle moving in elèctric field the equation takes the simple form 

dw dw „ d . . 

m-^ + Fv sx m = ( >- (87) 

Now we find thetomographic symbols of the integrals of motion. The kernels of the integrals of motion 
in position representation read 



x + 



2m 



iht 



5(x — x) ó'(x — x), 

m 



xo(x, x) - 

Po(x, x) = —Ft5(x — x 1 ) + ih5'(x — x'). 
Using these kernels one can find Weyl symbols of the integrals of motion 

W$ {q,p) = /^o (? + |,g-|) e-r^ u du = q-^t + 
Po(q+^,q- 2) e ~ ÏPU du = p- Ft. 



2m 



59) 
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The tomographic symbols of the integrals of motion are related the their Weyl symbols by the same 
formulas that the tomogram of a quantum state is related to the Wigner function (h = 1), i.e., 



W~ xo ( q ,p)e~^~^)^^ 

(2tt) 

<9 V 1 <9 t ( d \ ^ d Ft 2 



dXj d» m\dXj ^ + 2^' Wï(X ' /i ' Z/) ' 



dk dq dp 



d V 1 d 



dX dv 



+ Ft }w ï {X,fi, v) 



(90) 



The tomographic symbol of identity operator in these equalities has the form 



wï(X,h,v) 



-ikX 



k 2 



dk 5(p,)5(v) 



(91) 



The evolution equation for the tomogram can be solved by means of the Green function (propagator) 
of this equation relating the initial and current tomograms by the integral transform 



w(X, /i, v,t) = J U(X, /i, v, X', //, v', t)w(X', fi', i/, 0) dX' dfx' dv' . 
The propagator satisfies the evolution equation (we reconstruct the Plank constant) 



dU dU i 
~dt~ fi ~dv~ + h 



V 



d \ 1 d ihv d 
-7T- + 



dX) d fi 2 dX 



V 



d \ 1 d ihv d 



dX ) dfi 2 dX 



and the initial condition 



n = 0, 



(92) 



(93) 



U(X, fi, v, X', fi', v', 0) = S(X - X')5{fi - n')5{u - v') 



(94) 



It can be shown |15| [L6| that there is the relation of the Green functions of Schròdinger and tomographic 
evolution equations (h = 1) 



ïl(X, /i, v, X', //, v', t) = J k 2 G (^a + Y' y > f J G * ( a - Y > z > f j ~ z ~ ku "> 



x exp < ik i X' — X + [ia — /i 



dk dy dz da. 



For the charge moving in the elèctric field the evolution equation for the propagator reads 

<9I1 dU dU 

-m-^ + Fu dx = °- 



It has solution 



IL(X, /i, v, X', a', v, t) = 6(X - X' - Ft)5 \a-u! - —v) S( 

\ m J 



v — v 



(95) 



(96) 



(97) 
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Figure 3: Tomogram 



It is interesting to calculate the tomographic probability distribution of stationary state of the charged 
particle. Evaluating the corresponding integral fl80|) we get the tomogram of the state 



w £ {X, /i, v) 



1 



2ttM 



ifi 2 iX . 
2hv V -T^ V > dy 



Using explicit expressions for ipg one can obtain the explicit expression for the tomogram 



(9? 



w £ (X, /i, v) 



A 2 h 



iit\fi\ 



$ e- 



aX h 2 a 2 v 2 



16n 2 Fh\fx\ 



4 jj? 
f 2m 



F 2 h 2 



1/3 í2mF\ 1/3 X _ í m 2 F 2 \ 2/3 



(99) 



Plot of the tomogram is shown in Figure (|j). One can see the reach structure of the tomogram. 
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5 Conclusion 



Tha main result of our work is the obtaining the tomogram of the stationary state of the charged particle 
moving in the constant elèctric field. We found closed connection with integrals of motion, Wigner function 
and tomograms. The obtained results give a possibility to find new relations for the Airy function. 
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